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Abstract 

We study the distribution of Mahler's measures of reciprocal polynomials with com- 
plex coefficients and bounded even degree. We discover that the distribution function 
associated to Mahler's measure restricted to monic reciprocal polynomials is a recip- 
rocal (or anti-reciprocal) Laurent polynomial on [l,C)o) and identically zero on [0,1). 
Moreover, the coefficients of this Laurent polynomial are rational numbers times a 
power of TT. We are led to this discovery by the computation of the Mellin transform of 
the distribution function. This Mellin transform is an even (or odd) rational function 
with poles at small integers and residues that are rational numbers times a power of 
TT. We also use this Mellin transform to show that the volume of the set of reciprocal 
polynomials with complex coefficients, bounded degree and Mahler's measure less than 
or equal to one is a rational number times a power of tt. 

2000 Mathematics Subject Classification: 33E20, 44A05 

1 Introduction 

The Mahler's measure of a polynomial f{x) e C[x] is given by the expression 
(1.1) M(/)=exp|£log|/(e2"*)|dt|. 

If f{x) has degree M and factors over C as f{x) = wm Yl^^iix — /?«), then by Jensen's 
formula, 

M 



fJ-if) = \wm\ Yl max{l, 



m—l 



It is readily apparent that Mahler's measure is a multiplicative function on C[x]. In this 
sense Mahler's measure forms a natural height function on C[a;]. In this paper we study 
the distribution of values of Mahler's measure restricted to the set of reciprocal polynomials 
with bounded even degree and complex coefficients. 
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f{x) is said to be reciprocal if it satisfies the condition 



x"/(i)=/(x). 

If /(.'/;) Ih reciprocal and f{x) = Y^m=o'*^'mX™' , then it is easily seen that Wm = WM-m for 
TO = 0, ... , M. The reciprocal condition also imposes a condition on the roots of f{x): if 
/(a) = 0, then /(a^^) =0. If M = 2N there exists a Laurent polynomial 

N 

(1.2) p^{x)=Vo + Y,^n{x'' + X-^) 

n=l 

such that f{x) = x^Pv{x). We call Pv{x) the reciprocal Laurent polynomial with coefficient 
vector V. The collection of reciprocal Laurent polynomials with complex coefficients forms 
a graded algebra. 

The integral defining Mahler's measure makes sense for reciprocal Laurent polynomials, 
and it is easily seen that /x(pv) — IJ-if)- It is convenient to work with reciprocal Laurent 
polynomials since they form an algebra (the set of reciprocal polynomials is not closed under 
addition). We define the reciprocal Mahler's measure to be the function /Xrec : C''^"'"^ M 
given by 



(1.3) Mrec(v) = /x(Pv) = exp 




vo + 2 Vn cos(27rni) 



dt 



If V = {vq, . . . , vl, 0, . . . , 0) with vl 0, then there exist ai, . . . , a2L not necessarily 
distinct nonzero complex roots ofpv{x). By reordering if necessary, we may assume aL+n = 
a~^, and we may write 

L 

x^Pw{x) =vl '[[{x-Oin){x-a~'^), 

n=l 

and from Jensen's formula we have 

L 

(1.4) Mrec(v) = I^^lI Y[ max{|a„|, \a~'^\}. 

n=l 

From this expression we sec for all v S C^+^ and k G C the reciprocal Mahler's measure is: 

(i) non- negative: /ircc(v) > 0, 

(ii) homogeneous: /iroc(fcv) = |fc|/Xroc(v), and 

(iii) positive-definite: /Ltrec(v) =0 if and only if v = 0. 

In addition fi^cc is continuous as originally proved by Mahler [3] . 

By properties (i), (ii) and continuity, we find that firec is a symmetric distance function 
in the sense of the geometry of numbers (see for instance the discussion in [1, chapter IV]). 
(j,rec satisfies all the properties of a metric except the triangle inequality. The 'unit ball' is 
thus not convex. Explicitly, 

Vn+1 = {v e C^+l : Mrec(v) < 1} 
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is a symmetric star body. By property (iii) this star body is bounded. We call Vjv+i the 
complex star body determined by the reciprocal Mahler's measure. One of the principal 
results presented here is the computation of the volume (Lebesgue measure) of Vat+i. 
We introduce the monic reciprocal Mahler's measure, z^rec : K, defined by 

^rec(^) — Mrec 

Thus Uj.ec{h) is the Mahler's measure of the monic reciprocal Laurent polynomial 
(1.5) pb{x) = {x"" + x-^) +bo+Y,bn (x" + . 

n=l 

We denote Lebesgue measure on Borel subsets of C-*^ by A2JV, and introduce the distribution 
function associated with the monic reciprocal Mahler's measure, /ijv(0 • [O)0o) ~^ [0, 00), 
given by 

/iiv(0=A2w{beC^:J^rcc(b)<e}. 

/ijv(0 encodes statistical information about the distribution of Mahler's measures of recip- 
rocal polynomials with complex coefficients and even degree bounded by 2A''. 

The distribution function /lAr(^) is increasing and continuous from the right. By equation 

1.4 i/rec(t)) > 1 for all b e C^, thus /lAr(^) is identically zero on [0, 1). In fact, /ijv(l) = 0. 
To see this, suppose b G with frec(b) = 1. Then, from equation 1.4, Ph{x) has all its 
roots on the unit circle. Thus, if a is a root of Ph{x) then so is a = a~^. We find that 
b e , and hence the set of b € such that frec(t>) = 1 has A2jv-iiieasure 0. Thus 
/ijv(l) = 0, and /iAr(0 is continuous at ^ = 1. 

Wc recall the definition of the McUin transform. Given a function g : [0, 00) — > R, the 
Mellin transform of g is the function of the complex variable s given by 

We will give an explicit formula for /iAr(5) by computing its Mellin transform. We note 
that, since h]y{S,) is identically zero on [0, 1] the integral defining hN{s) can be written with 
domain of integration [1, 00). 

The integral defining hpfis) converges in the half plane 5i(s) > N. To see this, we use 
the following consequence of Jensen's inequality: 

M/)<ll/(^)l|2, 

where |l/(a;)||2 is the Euclidean norm of the coefiicient vector of f{x). Thus from equation 

1.5 we have 

z^rec(b) < (2 + |6o|' + 2|fei|2 + . . . + 2|6jv-ir)'/' < V^(l + l^or + . . . + \bN-i\Y^^ 
and hence 

{b e : z.,ec(b) < a C |b e : (1 + \bo\' + ... + l^jv-iH'/" < ^| ■ 
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The latter sot is a 'slice' of a solid sphere of dimension 2N + 1, and is thus a solid sphere of 
dimension 2N with radius less than ^/v^. Thus there exists a constant C such that 

/ijv(0 = A2jv{b e : i/,ec(b) < < Ce"" ■ 

It follows that 



hN{s)= J^^ r''hN{Oj<C ^' 



2N-2s ^ 



The latter integral converges if 3fJ(s) > A'', and hence hN{s) is defined in the half plane 
5R(s) > N. 

We follow the method introduced by Chern and Vaaler in [2] to express the volume of 
Vat+i in terms of the Mellin transform of ft.Ar(^). 

Theorem 1.1. For each positive integer N 

A2JV+2(Vjv+l) = 2TThN{N + 1). 

Proof. The volume of Vjv+i is given by 

(1.6) A2iV+2(ViV+l) = £ A2iV |b e : ^rec ^ ) < l| d\2{z). 

By the homogeneity of /Xrec we see 

A2Jv|beC^:/X,ec(^ ^ ) <l| = A2Ar jzC G : Mrec ) < 1 

= |^rA2;v|ceC^:Mrec( ^ ) < ^ 



JV 



and thus the integral in equation 1.6 can be written as 

JjzrhN (^^) dX2{z) = 27r j\^''+'hN(^l^ dr. 

The domain of integration in the latter integral is [0, 1) since hN{l/r) is identically zero on 
[1, oo). By the change of variables r = 1/^ we find 

/oo 
^_2(iv+i)-i/^^(^^ = 27rhN{N + 1). 

□ 

If we regard the integral defining /ijv(s) as a Lebesgue-Stieltjes integral, we may use 
integration by parts to write 



/lAr(s) = - 



2s 



1 

2s 
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Since /lAr(l) = and /lAf(C) is dominated by C^^^ , the first term vanishes when 5i(s) > A''. 
After a change of variables, we can write 

hN{s) = ^ ( I/rec(b)"^MA2Jv(b). 

The latter integral is interesting enough to name: 

Hn{s) = [ i/,ec(b)-2^dA2jv(b). 

The bulk of this paper is committed to the discovery that Hn{s) analytically continues 
to a rational function of s. 

Theorem 1.2. For each positive integer N, the function Hn{s) extends by analytic con- 
tinuation to an (even or odd) rational function. In particular, 



g2 _ 
n— 1 



Corollary 1.3. For each positive integer N, 

X2n+2{Vn+i) = (2iV + l)! ■ 

Proof. This follows immediately from Theorem 1.1 and Theorem 1.2. □ 

Corollary 1.4. For each positive integer N, /lAr(^) is a reciprocal or anti-reciprocal Laurent 
polynomial on the domain [l,oo) and identically zero on [0, 1). Explicity, if^>l then 



N 



n—l ^ ^ ^ ' 

Proof, /ijv(s) = Hn{s)/2s is a rational function whose denominator is a product of distinct 
linear factors of the form s — n. We use partial fraction decomposition to write 



We compute p{n): 



and so 



(s - n)/ijv(s) = JT 2^""^ 2 

s + n ^\ s"' — m"' 

m=l 



n-1 I ^ I 



Pin) = -^2--wn7;;^ n 



= 1 ^ ' m=n+l ^ > 



^ — vn?) 
n\ \( (-l)^-"(2n)! 



^ (n - l)!(2n -\)\)\{N^ n)\{N - n)\ 
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It is clear that p{—n) = (— l)^p(n), and so 
A quick calculation shows, for s > n, 

And so, by the uniqueness of the Mellin transform we find 

N 

(1.8) h^io = Yl 2pW(r'" + (-1) V") 

n=l 

for ^ € (1,00). The lemma follows by substituting equation 1.7 into equation 1.8. □ 

We outline the proof of Theorem 1.2. Given a G (C \ {0})^, wc can create the unique 
monic reciprocal Laurent polynomial Pa.{x) having ai, . . . , un, , . . . , as roots. We 
will use the change of variables a 1— > a to write Hpf{s) as an integral over root vectors of 
reciprocal Laurent polynomials, as opposed to coefficient vectors. This change of variables 
is useful, since by equation 1.4, v^cd^) is a simple product in the roots of 'Pa.ix) {i.e. in the 
coordinates of a). Analysis of the Jacobian of this change of variables will allow us to write 
Hn{s) as the determinant of an x TV matrix, the entries of which are McUin transforms 
which evaluate to rational functions of s. Theorem 1.2 will follow from the evaluation of 
the determinant of this matrix. 

Before proccding to the proof of Theorem 1.2, we present /ircc and "^^A^+i from another 
perspective. Given the positive integer Af, we define the Mahler's measure function to be 
fj, : C''^"'"^ R where /x(u) is the Mahler's measure of the polynomial with coefficient vector 
u. As was shown in [2], /i is non-negative, homogeneous, positive-definite and continuous. 
Thus fi is a symmetric distance function and the set 

UM+i = {ueC^+i:/x(u)<l} 

is a bounded symmetric starbody. Let M = 2N and consider the linear map A : C^+^ 
C2^+i defined by A(v) = {vo,vi, . . . ,vn-i,vn,vn-u . . . ,vi,vof . We define = A(C^+i) 
to be the subspace of reciprocal coefficient vectors. By equations 1.1, 1.2 and 1.3 we find 
Mrec(v) = ^(A(v)). Thus, thc starbody formed by the intersection 'U2Ar+i and V is related 
to the reciprocal starbody. Specifically, 

Vjv+i = A-i(T/nU2jv+i). 

Every bounded symmetric starbody uniquely determines a symmetric distance function [1, 
Chapter IV. 2 Theorem 1]. Thus, armed with and A, we could 'discover' /Urec- Equation 1.4 
can be recovered from the symmetry in the definition of A, so we would lose no information 
if we were to define fij-^c in this manner. 

The volume of Um+i, as well as the subspace volume of the starbody formed by inter- 
secting Um+i with the subspace of real coefficient vectors was investigated in [2] . Thus the 
computation of the volume of Vjv-i-i yields subspace volume information of another 'slice' 

of U2N+1. 
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2 A change of variables 

Let = C \ {0}, and define the map £jv : (C 

N 



C by £jv (") = a, where 



x^Pe^ix) = ]J(x + a„)(x + a„^). 



Thus the nth coordinate function of £;v is given by e„(Q;i, . . . , a^, a^^, ■ ■ ■ , a'^^), where e„ 
is the nth elementary symmetric function in 2N variables. Let : be the 

function whose nth coordinate function is e„, the nth elementary symmetric function in N 
variables. That is, given (3 e C-'^, if b = En {(3) then 

N N-1 
n—1 n=0 

It is well known that the Jacobian of -E'Ar(/3) is given by |y(/3)p, where 









/ 1 


1 


•• 1 \ 








Pi 


P2 . 


.. Pn 


(2.1) 


y(/3) = n (/^" 

l<m<n<N 


- Pm) = det 


P\ 


PI . 


.. P'n 
■■ P^-'J 



is the Vandermonde determinant. We will relate the Jacobian of £jv to the Jacobian of E^. 
Lemma 2.1. For each positive integer N, the Jacobian of En{oi) is given by 



V [ai + —,..., aN + — 

ai ajv 



n 

n 

n=l 



1 



Proof. By definition en{xi, . . . , xn , x[, . . . , x'j^^) is composed of all monomials of degree n in 



the variables Xi, 



.XN-X^, 



^n- 



If we impose the relation .x„i.x'„ = 1 for m = 1, 



.N. 



then Sn{xi, . . . ,xn,x'i, . . . ,a;^) is no longer homogeneous. In this situation it is easy to see 
that the monomials of degree n oi en{xi, . . . , xn , x[, . . . , x'^) consist of monomials which do 
not contain both Xm and x!^ for m = 1, . . . , A''. Hence, 

en{xi, ■ ■ ■ , xn,x'i, . . . , x'j^) = en{xi + x'l, . . . , xn + x'j^) + (monomials of degree < n). 

In general e„(xi, . . . , .tjv, x'l, . . . , x'j^) contains monomials of degree n — 2Ai formed from 
monomials which contain Xm and x'^ where m runs over a subset of 1, . . . , iV of cardinality 
M. By counting the number of times each monomial of degree n — 2M appears we arrive 
at the following identity. 



ai, . . .,aN, ■ 



ai 



1 

OIN 



ai 



+ 



+ 



N - 
N- 

[N/2] 



1 

- 2 



= E 



M=0 



- n 
1 

- n 
2 

N -n-2M 
M 



1 \ 






. . . , Q!iv H 












( 1 






en-2 Oil ^ , . 






V Oil 






( 1 






e„_4 ai H , . 


. ,aAr H 




V "1 







+ ... 



-2M ai 



1 

ai '' 



. . . , Q!Ar + 



OLN , 
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where [N/2] is the integer part of N/2. 
Thus 



/I ... 0\ 



£iv(a) 



* 1 



where 



0L\ 



■ 1/ 



and * represents entries which are not necessarily 0. The Jacobian of £^{(3) = \V{/3)\'^, 
and thus by the chain rule we arrive at the formula for the Jacobian of Sat (a) given in the 
statement of the lemma. □ 

The Jacobian of Eat (a) is nonzero for A2Ar-almost all points of (C^)^, and there are 
2^N\ preimages for A2jv-almost all a e C^. Employing the change of variables formula, we 
find 



Hn{s) 



/ fc'rec(a) ^''(iA2Ar(a) 

' N 

n max{|a„|, \a~^\} 



2^N\ 



(Cx)' 



(2.2) 



V{ai + —, 

ai 



dX2N{oc). 



The latter integral admittedly looks formidable, however this change of variables is beneficial 
since free is a simple product. 



3 H]sf{s) is a determinant 

We first prove a short technical lemma concerning determinants. 

Lemma 3.1. Let N be a positive integer. If I = I{j, k) is an N x N matrix and is the 
Nth symmetric group, then 

1 ^ 
(3.1) det(7) = ^ E E sg'iW 'Sr^(^) U ^(^W' ^W)" 

rGSiv ceSiv n=l 

Proof. 

N N 

l[l{T{n),a{n)) = I{n,a o T-'{n)). 

n=l n=l 

Thus we can write (3.1) as: 

1 _ ^ _ 1 ^ 

— EE sgn((Tor ^)Y[l{n,aoT \n)) = E E sgn{a)Y[l{n,a{n)) 

tGSjv <7&Sn n=l tGSjv <7&Sn n=l 

N 

= E ^^^^^'> n Iin,Hn)), 

(^ESn n=l 
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which is the famiUar formula for det(/). 



□ 



Using equation 2.1, we expand the Vandermonde determinant as a sum over the sym- 
metric group to find 



1 1 
V [ai-\ ,...,a„-\ 



^ f 1 



which we rewrite as 



N 





2 \<^(")-i / 




(an H 


h — I a„ H 


^ an) 




an J V 





cr(n)-l 



T(n)-1 



ctESn tESk n=l 

Substituting this expression into equation 2.2, we can write iljv(s) as 

' N / , , X 2N 

Jl max{|a„|,|a-^|}" 



1 



2^7V! 



N 



1 — ' 




a„ H 









T(n)-lN 



dA2jv(Q:). 



\cr£5N tESk n=l 

Exchanging the sums and the integral, and consolidating the products, we find 



Hn{s)= Y sgn(o-)sgn(T)-j^ / i f[max{|a„|,|a„^ 



a„ 



Y 

(n)-l 



-2s 



dX2N{OL). 



By an application of Fubini's Theorem we find 

(3.2) Hr,{s) = ^ E E ^^^^^) ^S'^M ft J('^W'^W) 



n=l 



where 3 {J, K) is given by 
1 



max 



{\al\a-'\] 



-2s I l\ _ l\ 1 

a \ a — — \ aH 

a I \ a I \ a 



J-i 



_ , l\^-'dA2(a) 



Applying Lemma 3.1 to equation 3.2 we find i?jv(s) is the determinant of the N x N matrix 
3 = 3{J,K). 



4 The entries of 3 are rational functions of s 

We shall view 3{J,K), not only as an entry in a matrix, but also a function of s. We note 
that \2(a)/\a\'^ is normalized Haar measure on C^. Thus 3{J,K;s) is invariant under the 
substitution a a~^, and we may write 



3{J,K;s) 



\a\ ^'^ I a - 



a — — ] an 

a J \ a I \ a 



j-i 



1 \ ^ ^ d\2{a) 
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where D is the open unit disk. By setting a = re'^ we may write 3{J,K;s) = h{J,K;r), 
where h{J,K;r) is given by 

for r £ [1,00), and identicaUy zero on [0, 1). 

By the change of variables 6 —6 we see that h{J,K;r) = h{K,J;r). We conclude 
that J is a symmetric matrix whose J, K entry is h{J, K; s). 

Lemma 4.1. 3{J,K;s) analytically continues to a rational function. Specifically 

n=l 

where 

' -j-u t J-l 



Cn{J) 



(';+„) — (.7+r._i) if n < J and n = J (mod 2) 

2 2 ^ i 

1 if n = J 

otherwise 



Proof. Without loss of generahty we assume K > J. There is a constant 6 such that 
h{J,K;r) < Gr^~^^ on [l,oo). Thus the integral defining 3{J,K;s) converges in the half 
plane 3?(s) > {J + K)/2. 

Writing (re*^ + and + ^)^~^ as sums with binomial coefficients we may 

rewrite 4.1 as 

(4.2) hiJ,K;r) = E E ; ') (""^ ') 

X J^^ + ^ - (e"'" + e""'") ) e(^-^-2(^-'=»^^ d0. 
The integral appearing in this expression can be readily evaluated: 

J^^ (r^ + 1 - (e^*^ + g(J-/f-20-fc))ie 

r 27r(r2 + ^) fc = j + (if-J)/2 

(4.3) = -27r fc=j + l + (ii'- J)/2 . 

[ -27r k= j -1 + {K - J)/2 

If J ^ K (mod 2) wc sec that h{J, K: r) (and hence 3(,/, if; s)) is identically zero. 

The conditions given in equation 4.3 allow us to eliminate one of the summations in 
equation 4.2. We use the facts that Q < k < K — 1 and < j < J — 1 together with the 
conditions in 4.3 to find conditions on j. Specifically, 

k = 3 + => < j < J - 1, 

fc = i + l + ^ 0<i<min|:^-2,J-l|, 
k = j-l + ^^^ max|^^^ + l,o| <j< J-l. 
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Since K > J, we can write 

J - K 



max 



1,0 



and 



J + K 



2, J - 1 



J-1 



where 6jk = 1 if J = K and is otherwise. From this information we may write h{J, K; r) 



as 



J-1 

3 



K- 1 

K-J 



+ 3 



J-1 

- E 



J-1 

3 



K-1 



„2J-4j 



J-1 

E 

J-l-Sj!. 

■ E 

i=o 



J-1 



ii'- 1 

K-J 



+ 3 



„2J-4j-4 



J-1 



K-1 



„2J-4i-4 



Using the convention that {^^^) = and = we may eliminate 5jk from the latter 

two sums. Reindexing each sum based on the powers of r which appear and simplifying the 
binomial coefficients, we find 



h{J, K-r) = 27r ^ 



J -1\(K~1 



J -1\ fK -1 



+ 1 \^ + 



„4i 



E 



+1 



J-1 



K - 1 



r> ^ / V O ^ 1 



„4; 



iC- 1 



Note that in the case that J is odd, these sums run over consecutive odd multiples of 1/2. 

Reindexing the first and third sum by ? — ^, we may combine the first and second sums, 
and the third and fourth sums. We may then write h{J, K; r) as 



K- 1 

+ 1-1 



Due to the symmetry in the summands we may reindex the sums using only positive indices. 
Let /o = if J and K are even, and Zq = 1/2 if J and K are odd, then 



h{J,K;r) 



2n 



K-1 



J+K 



- 1 



K-1 

J+K 
2 



+ 27r^ 



J-1 



J-1 

i+l-1 



K - 1 

f +/ 



K-1 
f +/-1 



We are now in position to compute /i( J, K; s). There is a correspondence between the co- 
efficients and powers of r which appear in h{J, K; r) and the poles and residues of h{J, K; s). 



As was demonstrated in the proof of Corollary 1.4, the Mellin transform of 



analyt- 



ically continues to the rational function s/{s — 41 ). Thus 3(J, iir;s) extends to a rational 
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function: 



3{J,K-s) 



2tt 



K-l 

J+K _ 1 
2 



+ 



J -I 

i+l 



K-l 

J+K 
2 



J- 1 



J2 



K - 1 



if - 1 
+ 1-1 



K 



72 ■ 



Reindexing this sum by setting n = 21 we find 

' K-l \ /K-l 

J-1 

J+n 
2 



'i{J,K;s) = 277 

+ 2.^ 



2 

J-1 



- J2 



J+n _ 1 
2 



K - 1 

K+n 
2 



K-l 

^-1 



where the sum is over n G {1, 3, . . . , J — 2} if J and K are odd. and over n € {2, 4, . . . , J — 2} 
if J and K are even. If J = if the leading coefficient is 1. If we set 



Cn{J) 



then we may write 



(i+l) - ( j±^_i)J if n<J and n=J (mod 2) 

if n = J 
otherwise 



N 



(4.4) 



3{J,K;s)=nJ2cn{J)cn{K)- 



2s 



It is easy to verify that this expression is symmetric in J and K, giving 3{J,K;s) = 
3{K,J;s) as expected. Additionally if J ^ if (mod 2), the expression in equation 4.4 

yields 3{ J, K; s) — 0. This proves the lemma. □ 

We identify J( J, if; s) with the rational function it extends to. When J and if are odd, 
J(J, if;s) has poles at ±1, ±3, . . . , ± min{ J, if}. When J and if are even, J(J, if;s) has 
poles at ±2, ±4, . . ., ± min{ J, if}. 3( J, if; s) has a zero of multiplicity one at 0. 

We are now in position to prove the first part of Theorem 1.2. iijv(s) is the determinant 
of J, and the entries of J extend to rational functions of s. Since the determinant is a 
polynomial in the entries of a matrix, Hj\[{s) itself extends to a rational function of ,s. In 
fact, since the determinant is a homogeneous polynomial in the entries of a matrix and the 
entries of 3 analytically continue to odd functions, Hn{s) analytically continues to an even 
rational function when N is even, and analytically continues to an odd rational function 
when N is odd. We also see Hn{s) has a zero of multiplicity N at 0. 



5 Hn{s) is a simple product 

In this section we express det(3) as a simple product. The structure of the poles and residues 
of 3{J, if; s) will allow us to find linear dependence relations on the rows of 3. 
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Let Bn be the N x N matrix whose J,K entry is the integer c„(J)c„(-ft'). Then by 
Lemma 4.1 we have the matrix equation 



N 



n=l 



2713 



Define e to be the row vector given by a; J = {cn{K))^^-^. It follows then that the 
Jth row vector of B„ is given by c„(J)wJ, and thus every row of B„ is a scalar multiple of 



We may find a nonzero vector tj) G such that = for 1 < n < A/' — 1. In fact, 

Bntp = for 1 < n < N — 1, leading us to the vector equation 



N 



n=l 



N2- 



We see that - B^-^tzw^ V' = 0) and so det ^3 - B^^rzw^ = ^- From the definition of 
Bn and Lemma 4.1, we find 



3-B 



2-KS 



/a(l,l) 3(1,2) 
J(2,l) J(2,2) 



a(i,7v) 

3(2, A^) 



V3(Ar,l) J(7V,2) ... 3{N,N)-^^J 
Taking determinants and exploiting the multilinearity of the determinant 



det 



/a(l,l) 3(1,2) ... 0{1,N)\ 
J(2,l) 3(2,2) ... J(2,Ar) 



V3(^,l) a(iV,2) ... 3{N,N)J 



det 



/3(1,1) 3(1,2) 
3(2,1) 3(2,2) 

V3(7V,1) 3(iV,2) 



\ 




The left hand side is Hn{s). By a simple induction argument we finally arrive at a simple 
product formulation of Hf^{s): 



N 



hn{s) = n 



2tts 
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